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I review recent progress made on heavy quark physics on the lattice. 



1. INTRODUCTION 

Precise knowledge on B meson decay properties 
plays an essential role in testing the unitarity of 
the CKM matrix The experimental status 

relevant to the unitarity test is now promising be- 
cause experimental uncertainties are already very 
small or expected to be reduced well in the near 
future ||. On the other hand, theoretical cal- 
culations relevant to the test are still not suffi- 
ciently accurate due to the non-perturbative ef- 
fect of QCD. Lattice QCD is an ideal tool to deal 
with this effect and should be able to reduce the 
current theoretical uncertainties Q]. 

Apart from the asymmetric e + e~ colliders, 
there are two important upcoming experiments 
for the unitarity test. The Tevatron Run II will 
produce a large number of and variety of b and c- 
hadrons, and their basic properties such as mass 
and lifetime will be precisely measured || . What 
is important for us is that once B®-B® mixing is 
observed the mass difference AMb s will be mea- 
sured with a few percent accuracy. The width dif- 
ference in the B s meson system is also expected 
to be measured precisely. Another exciting exper- 
iment is the CLEO-c project ||. There, charmed 
quarkonia, hybrid states and glueballs will be ob- 
served. In particular the leptonic and the semi- 
leptonic decays of D mesons are expected to be 
measured with a few percent level. Advancing 
lattice QCD calculations with a view to combin- 
ing them with these precise experimental results 
is an urgent task in front of us. 

In this review, I will focus on new and updated 
calculations of hadron matrix elements. The cur- 
rent status and progress made in spectrum cal- 
culations and formulations are not covered, al- 
though they are very important and interesting. 



In particular, I will mainly discuss the lattice de- 
termination of the B°-B a mixing amplitude and 
how to put a strong constraint on the poorly 
known CKM element \V t d\- The other important 
quantities such as form factors of semi-leptonic 
decays will be mentioned briefly. Recently sev- 
eral suggestions to improve the limits on accuracy 
of present lattice calculations have been made 
in methodology. I will briefly introduce some of 
them before summarizing my point of view on the 
current status. 

2. B°-B° mixing 

2.1. General remarks 

Within the Standard Model, the mass differ- 
ence of neutral B meson system is given by 

AM Bq = (known factor) x \V t * b V tq \ 2 

(B° q \h„(l-75)q Ml -75)9 
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M Ba 
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where q=d or s. AM Bd has already been mea- 
sured accurately AM Bd = 0.503 ± 0.006 ps" 1 
, while for AM Bs only a lower bound is known 
(AMb, > 14.4 ps" 1 at 95 % CL §), but it is 
expected to be measured precisely at the Teva- 
tron Run II. The hadron matrix element in (|l]) is 
usually parameterized as 



(B° q \ h„{l-l 5 )qh^ 



-Jl a B B U b )M_ 



where the decay constant f Bq 
the following matrix element, 

\B°„) 



-7 5 )<Z \B° q ) 

(2) 

is defined through 



M 1 " 75)9 \B v q ) = if BqPll . 



(3) 



Since the four-quark operator receives renormal- 
ization, B B depends on renormalization scheme 



year 


f Bd [MeV] 


B B d 


'97 Onogi 


163(12) 




'98 Draper 


165(20) 


1.32(6)(12) 


'99 Hashimoto 


170(20) 


1.23(23) 


'00 Bernard 


175(20) 


1.30(12) 


'01 Ryan 


173(23) 


1.30(12) 



Table 1 

The quenched f Bd and B Bd quoted in the recent 
reviews. 



and scale. In the literature, the renormaliza- 
tion scale independent B parameter is often used, 
which is defined by 
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where the full expression of J n , is found in ||. 
Throughout this review, I quote B B rather than 
the scale dependent one B Bq (aO ■ 

In the following, I summarize the current status 
of the quenched and unquenched calculations of 
f Bq and B Bq . 

2.2. Quenched f Bq and BB q 

According to the recent Lattice conference re- 
views P, p0| , pTp^Jl3] | , the quenched results for f Bd 
and B Bd have been stable over several years as 
shown in Table @. This year JLQCD @ up- 
dated their quenched calculation, which is per- 
formed at /3=6.0 using the non-perturbatively 
0(a) improved Wilson light quark and NRQCD 
for heavy quarks. They find f Bd = 158(3)(10) 
MeV, f B Jf Bd = 1.16(1)0, B Bd = 1.29(5)(7) 
and B B jB Bd = 1.020(24) (15) (±£), where the 
first error is statistical, the second systematic and 
the third comes from the uncertainty of m s . 
Figure [j] summarizes recent results for 
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quenched f Bd y,^,|lj,y,^,£2],£3| 
including the JLQCD's updated result [14]. The 
results are grouped into three categories, "FNAL" 
H), "NRQCD" (27) and "NPSW", depending on 
the formulation for heavy quarkQ We find an 
agreement within about 10% among the results 
obtained with the three different actions. Since 
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see 11 



brief summary on the actions used, for example, 
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Figure 1. Summary plot of quenched f Bd . Shaded 
band represents the current average given by 
Ryan CI. 



the new result by JLQCD is consistent with the 
previous average given by Ryan |l3| (shaded band 
in the figure), I keep her values as the summary 
as of this conference. 



f (N,=0) 
J B d 
f (N f =0) 
J B, 



,(N f = 
JDd 



173(23)MeV, 
200(20)MeV, 
203(14)MeV, 
230(14)MeV. 



(5) 



These values are helpful in calibrating new for- 
mulation of heavy quarks. 

As for B Bq , only a couple of calculations are 
available. Figure || shows the l/Mp d dependence 
of B Bd obtained by APE [§8j with the relativistic 
heavy quark (NPSW) and by JLQCD Q with 
NRQCD for heavy quark. UKQCD has also used 
the relativistic heavy quark p5|, but since their 
result agrees with the APE result very well it is 
suppressed in the figure for simplicity. I also plot 
the result calculated in the static limit (HQET) 
by Gimenez and Martinelli (29)0. 

The 1/Mp d dependence of B Bd seems incon- 
sistent between two formulations of heavy quark, 

2 Later their result was corrected in |jCJ|. 
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Figure 3. Summary of quenched Bb- Shaded 
band represents the current average. 



Figure 2. Comparison of 1/Mp d dependence of 

B B d - 



but the results agree within error at the physi- 
cal point l/M Bd =0.19 GeV -1 . A possible reason 
for this disagreement is systematic errors, which 
are different for different formulations; the rela- 
tivistic action has sizable 0(a 2 m 2 ) errors and the 
NRQCD approach contains large 0(a 2 s ) errors in 
the matching. In |3l), the results with the rel- 
ativistic approach is combined with the one ob- 
tained in the static limit and interpolated to the 
physical Bd meson mass (filled triangle up). It 
is interesting that the result with the combined 
analysis is slightly lower than that with extrapo- 
lation (filled triangle down) and the consistency 
with the NRQCD becomes better. 

SPQcdR has started a new calculation at 
/3=6.45 H with the RI-MOM scheme @ for 
renormalization. The strategy is similar to what 
APE adopted, namely relying on the extrapola- 
tion in heavy quark mass from around c quark re- 
gion to the the physical b quark mass. However, 
due to the use of high (3 value, the distance to the 
physical B meson mass becomes smaller. The mo- 
tivation of this work is to test if the 0(a 2 ) error 
is under control. While the numerical results are 
still very preliminary, no significant change has 
been observed compared to the APE results at 
/3=6.2 [^8|. This suggests that the discretization 
error, which is most serious in this extrapolation 
method, is under control. 

Figure [?] summarizes the quenched calculations 



of Bu d - At this moment, I quote the following 
conservative value as a current estimate, 



^ Nf " 0) = 1.33(12), 



(0) 



which is indicated in the figure by shade. 



2.3. Unquenched /n 

Several large scale unquenched calculations of 
Jb were carried out a few years ago, and the re- 
cent reviews reported 10-20 % increase compared 
to the quenched ones |i^ , |lT| , ^2|Jl^] . I wish to dis- 
cuss that this conclusion needs a reexamination 
because of possible effects of chiral logarithms ex- 
pected from chiral perturbation theory 

Incorporating heavy quarks into chiral pertur- 
bation theory via heavy quark symmetry was 
first proposed by Wise Q and Burdman and 
Donoghuc [f35f . Since then several authors have 
developed this idea jjjB. Recently Booth (3?) 
and Sharpe and Zhang |38| extended the idea 
to (partially) quenched chiral perturbation the- 
ory ((P)QChPT) to make it applicable to actual 
lattice calculations. 

In the Nf~2 unquenched case, where u and 
d quarks are treated dynamically, PQChPT pre- 
dicts the non-analytic term in fs d to be p8| 



r(l + 3 3 2 ) 



(4tt/)- 



ln^/A 2 ), 



(7) 



where A is the cutoff of the theory and /= /fl- 
at the order considered. While /^ has a similar 
expression for its non-analytic term, in this case 
the B*Btt coupling <?, which is only poorly con- 
strained in <7=0.3~0.6 [Bh, appears because the 
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Figure 4. Chiral extrapolation of $/ B given by 
]39fl. Both axis are normalized by Sommer's scale 

mass difference in £? and £?* mesons is smaller 
than the pion excitation. For fs B , the expression 
of the non-analytic term takes a different form. In 
the Nf—2 partially quenched case, where u and 
d quarks are treated dynamically as before and s 
appears only as a valence quark, the prediction 
for the non-analytic term becomes 



(l + 3<? 2 ) 



(4tt/) ; 



Hm%/A 2 



(8) 



where m 2 g = (m 2 s + m 2 T )/2, We observe that 
while for the chiral behavior of fs d (0) indicates 
significant slope as m 2 vanishes, suggests a 
milder behavior for fg a because rnr s ^ becomes a 
constant when m 2 s 3> m 2 . Lattice data are sup- 
posed to show these logarithmic behaviors. 

JLQCD has accumulated twice as many statis- 
tics as they had last year in Nf=2 QCD, and 
updated their analysis of the decay constant [39] . 
The simulation is performed at /3=5.2 (a ^0.09 
fm) with the NRQCD action for heavy quark 
and the non-perturbatively 0(a)-improved Wil- 
son quark action for both dynamical and valence 
light quarks, which range from 0.5 m s to 1.5 m s . 
The chiral extrapolation of $ f B ~fB q yjMs t is 
shown in Fig. |[ They attempt to fit 5>/ B and 
<L>/ Bs with (0) and (|sj) , respectively, as well as with 
a conventional quadratic fit. They find that fs d 



depends on the fit form and the B*Bir coupling 
g significantly, while fs s does not, in agreement 
with the observation made above. This means 
that /s 3 //s d has large uncertainties as well. 
From this analysis, they give /s j> //B d =1.24-1.38 
as a preliminary result. 

Recently, Kronfeld and Ryan jidj] pointed out 
the similar enhancement of the ratio by taking 
into account the predicted chiral log and recon- 
sidering the chiral extrapolation. Their estimate 
is W/ B<J =1.32(8) and B B ,/B Bd = lM{2). Both 
analyses suggest that the ratio can be signifi- 
cantly larger than the previous world average 
/b s //b,=1.16(5) @. 

Here let me summarize the Nf=2 unquenched 
calculations of fa. At present, consistency be- 
tween lattice data and ChPT is not clear because 
the value of g is still unknown. In addition, the 
0(1/M) contribution to the log term is not known 
well [^2j . In this review, the central value of de- 
cay constants is taken from the previous ones |Ts)| , 
which can be considered as those obtained with 
a quadratic fit in the chiral extrapolation. The 
uncertainty for fs d associated with the chiral log- 
arithm is estimated following JLQCD while it is 
neglected for fg, ■ According to their analysis, 

(f^7-fEf T )/f^T = 0.17, (9) 
where /^" a is obtained by quadratic form (as 
usual) and /s' lPT includes the effect of chiral log- 
arithm. In obtaining /g hPT , they take the upper 
bound <7=0.6, for which the effect becomes max- 
imum. Taking this uncertainty into account, I 
quote 

/^ f=2) = 198(30) (±j°) MeV, (10) 
j(N t =2) = 230(30) MeV, (11) 



(N f =2) , , (N f =2) 



1.16(5)(3 



+24\ 



(12) 



as my conservative estimates. 

The physical prediction of fs d s requires three- 
flavor dynamical simulations. MILC has started 
such an attempt |Q. They employ a highly 
improved gauge and staggered quark actions for 
gauge configurations and the tadpole-improved 
clover action for valence quarks, where FNAL for- 
malism is applied for heavy quarks. The simula- 
tion is performed for two lattice spacings a ~ 0.13 
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Figure 5. Chiral extrapolation of af Qq in valence Figure 6. Chiral extrapolation of f Bd in dynami- 

quark mass given by ||. ca l quar k mass g i ven by E3. 



and 0.09 fen. To keep the lattice spacing constant, 
(3 is adjusted as dynamical quark mass is varied. 
Since Za is not available yet, they focus on ratios. 
Chiral extrapolations are made in two steps. First 
the valence light quark is extrapolated (interpo- 
lated) to the physical up/down (strange) quark 
mass at each dynamical quark mass. Then dy- 
namical quark mass is extrapolated to up/down. 
The third flavor of dynamical quark is tuned to 
m 8 in advance. 

Figure ^| shows the chiral extrapolation of afQ q 
in valence quark mass. They attempt to fit the 
data by two types of fit forms, linear (solid line) 
and quadratic (dotted line) functions, and find 
no significant deviation between them. It is in- 
teresting to see whether this linear behavior is 
consistent with Nf=3 PQChPT. But such a test 
is complicated because of flavor breaking effects 
in pion masses for the staggered quark action Q . 

The chiral extrapolations of fs d and /s //s d in 
dynamical quark mass are shown in Figs. j|and 0, 
respectively. In || where the scale is set by m p , a 
negative slope is seen. But once one sets the scale 
by r\ f45|] , the slope becomes unclear, namely fg d 
appears to be independent on md yn - In 0, they 
find that the ratio is independent both of Nf and 
TO d yn, and quote /s s //s d = 1.18(1)(^) (prelim- 
inary) . It is very interesting to check if the whole 
above observations are consistent with ChPT. 



2.4. unquenched Bs q 

This year no new or updated calculation on the 
unquenched i?-parameters was presented. Last 
year JLQCD reported that the fit form depen- 
dence is relatively mild for Bs q in their Nf = 2 
unquenched simulation |^6|. This is consistent 
with the PQChPT prediction that the chiral log- 
arithm term for BB d is given by Q as 

1 777 ^ 

--(l-3 5 2 )^^lnK/A 2 ), (13) 

and hence suppressed compared with the case 
of fs d in (0), and that chiral logarithm is al- 
together absent for Bb s - JLQCD also reported 
that quenching effects are small from a compari- 
son with the quenched Bs q ■ 

The ratio BB S /Bs d is also stable against the 
number of dynamical flavor, the heavy quark ac- 
tion and the fit form in the chiral extrapolation, 
and the recent results are compared in Fig. |^, 
where only the data shown in top is unquenched 
results. 

Since at present there is only one unquenched 
calculation by JLQCD, I take their values as cur- 
rent estimates, and quote 

< f=2) = 1.34(13), (14) 




dynamical (m^ / m p ) 



Figure 7. Chiral extrapolation of f B J f Bd in dy- 
namical quark mass given by [}43| . 



B 



(N f =2) 



D 



(N f =2) 
-B rf 



= 0.999(24). 



(15) 



It should be noted that the ratio is already very 
precise. 

2.5. SU(3) flavor breaking ratio 

The Tevatron Run II experiment is expected to 
yield a precise value for the mass difference of B s 
meson system, AM S , in the near future S . Once 
AM S is determined precisely, assuming the three 
generation unitarity (|Vf S | |^c&|) we can put a 
strong constraint on the poorly known CKM ma- 
trix element \Vtd\ through the SU(3) flavor break- 
ing ratio of the hadron matrix element, 



AM F 



\V ts \ 2 Mb. 



\Vtd\ 2 M Bd 



(16) 



where £ = (/ B . y 'B B ,)/(f Bd y/ 'B Bd ). 

Until recently it has been expected that £ can 
be easily determined within a few percent accu- 
racy by lattice QCD. However, the presence of 
chiral logarithm make s it more difficult. As we 
discussed in Sec. 



2.3 



fB s /fB d can still have a 
large systematic uncertainty while the ratio of 
B Bq does not. Following the values given in (|l2| ) 
and (151), I quote as my personal estimate, 



£ 



(N f =2) _ 



= 1.16(6)(^ 4 ) 



(17) 
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Figure 8. Comparison of quenched and un- 
quenched B B jB Bd . Only the top one is ob- 
tained from an unquenched simulation. The cen- 
tral value I quote is indicated by shade. 



Kronfcld and Ryan also gave a similar value 
(£=1.32(10)) in their analysis 

I now discuss that a precise determination of 
| Vtd | is still possible if one considers the Grinstein 
ratio of decay constants defined by 

Rl = (fBjfB d )/(f D JfD d ). (18) 

Rewriting (|l^) in terms of R±, one obtains 
\V ts \ 2 (M B \ 



( AM Bt 
V AM B . 



Runll 



\V t d\ 2 \M Bd J known 



4* 



B f 



B B d 



.(19) 



Lattii 



CLEO-c 

As indicated in the subscript, the ratio /d,/ 1 fo d 
is expected to be measured in CLEO-c, and 
B B jB Bd can be determined by lattice QCD pre- 
cisely. Therefore R\ is the only remaining quan- 
tity to be fixed. 

We expect this task is achievable in lattice 
QCD since variations under chiral and heavy 
quark expansions will both cancel out in the ratio 
to a large extent, leaving only a small deviation 
from unity. 

This year JLQCD presented their preliminary 
result of R\ Jlsfl , The simulation is carried out 
on Nf=2 dynamical configurations at /3=5.2. In 
order to handle the c quark, FNAL formalism is 
applied for heavy quar ks. Figure p| shows the chi- 
ral behavior of R x x y/(M Bs /M Bd ) /(M Ds /M Dd ) . 
Several fit forms are attempted as shown in the 
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Figure 9. Chiral extrapolation of R\ x 
y/(M B jM Bd )/(M Da /M Dd ) given by @. 



figure. They found that the result is relatively 
insensitive to the fit form and its dependence is, 
at most, about 1%. Their preliminary result is 
Ri = 1.018(06)(10), where the first error is sta- 
tistical and the second represents the systematic 
uncertainty. 

In Fig. I gather the Grinstein ratio obtained 
by previous calculations. Only FNAL'98 |L5j and 
JLQCD'02(clover) §§ explicitly calculated R x . 
For other results, I only show their central values. 
It can be seen that this ratio is very stable and 
there is only a negligible fluctuation. At present, 
I quote 



(N f =2) 



1.02(2). 



(20) 



From a comparison with quenched and un- 
quenched (JV/=2) results, it is unlikely that the 
central value and its accuracy change drastically 
when one goes to Nf=3. 

As we have seen, once AM Bs , /d s a- n( i Id a are 
measured in the forthcoming experiments, we can 
test the unitarity at a few percent level through 
|Vtd|. To gain more confidence, the current val- 
ues of i?i and B B JB Bd should be confirmed by 
several groups. 

3. Form factors 

3.1. heavy to heavy transition 

In order to test the CKM unitarity, a precise 
determination of \V c b\ is indispensable, as it sets 
the normalization of the parameters (p,rj). To 
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Figure 10. Summary of the Grinstein ratio. 



this end, B — > D*lvi is most promising. The 
differential decay rate is given by 

— = (known factor) \V cb \ 2 \F B ^ D *{u)\ 2 , (21) 
aw 

where u> = v' D » ■ v B . At zero recoil (lu = 1), 
T B ^D'{1) = h Al {\), where ^(1) is given by 



(D*(v)\c W b\B(v)) = W2m B 2m D *e' li h Al (l). (22) 

Last year Fermilab group demonstrated a precise 
determination of T B —,d* (1) by a full use of heavy 
quark symmetry and obtained T b ^d* (1) = 
0.913lHm±o:o302 in the quenched approxima- 
tion p9[ . In my opinion the lattice method for 
calculation of this decay mode has been estab- 
lished. The extension to unquenched calculations 
remains to be done. It is worth noting that chiral 
logarithm terms would not affect this accuracy 
because they are suppressed by (1/M) 2 . 

At this conference a new calculation of A& — > 
K c lv was presented by Tamhankar |5C|] . This 
process is experimentally challenging, but gives 
an independent determination of \V c b\- They use 
the quenched 20 3 x64 lattice with cr x =1.32 GeV. 
The actions used are 0(a 2 , a s a 2 ) improved gauge 
and tadpole improved clover and FNAL approach 
is taken for b and c quarks. The chiral limit is not 
taken at present. Instead, they investigate the 
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Figure 11. The current status of quenched calcu- 
lations of the B — ► nlv form factors. 



initial heavy quark mass (mb) dependence with 
fixed m c and mi; g ht. No clear mass dependence 
has been observed in their preliminary analysis. 

3.2. heavy to light transition 

Once | | / 1 V^t, | and angle </>i(/3) are deter- 
mined precisely enough, the location of the apex 
is essentially fixed in the p-rj plane. Then the next 
step is to test the consistency of the CKM mech- 
anism, for example, by measuring other CKM el- 
ements such as \Vub\- B —> nlv is one of the sim- 
plest choice for such a purpose. The definition of 
form factors, / + and / , is given by 



(7r(p') I q^b | B(p)) 



+ 



(23) 



where q = p — p' . Figure |ll| [|5lJ shows the 
quenched calculations made by four major groups 
|5^,^,|^j5|] . Data from different groups agree 
with each other, but only within a large uncer- 
tainty (*~20%). Unquenched simulations remain 
to be done. Moreover the issue of chiral loga- 
rithm could be a significant source of uncertainty 
in this case. More studies are needed before fi- 
nalizing the form factor calculation. 

This year SPQcdR made two contributions on 



the heavy to light vector meson semi-leptonic de- 
cays, B -► plu (D -► K*lv ) @ and B -> K*-y 
pq . Both are performed at (3=6.2 and 6.45 with 
non-perturbatively 0(a)-improved Wilson action 
for light and heavy quarks, and the currents are 
also improved non-perturbatively. The form fac- 
tors of these decays obtained around c quark mass 
is extrapolated to the physical b quark mass by 
using the HQET scaling laws for a fixed value of 
v ■ p' or g 2 =0. According to their preliminary re- 
sults, no scaling violation has been seen for both 
matrix elements, but the statistical uncertainty is 
still significant. At present, chiral extrapolations 
are carried out by assuming the form factors to be 
a linear function of light vector meson mass. But 
in the future a guiding principle will be necessary 
to make a extrapolation more reliable. 

3.3. Determination of g^BTr coupling 

The determination of the coupling constant g 
defined in ChPT, or equivalently gp*p^ given by 

(P(p) n(q) | P*(p',X)) 

= -gp-p*tf)Q'^(p')P*)W-P-Q), ( 24 ) 
gp*Pn = lim o gp*Prr(q 2 )- (25) 

q 2 — 

is now extremely important because this coupling 
plays a crucial role in the chiral extrapolation of 
decay constant and in the normalization of one of 
the P — ► 7r form factor f + (q 2 ) at q 2 =q 2 aax - gp' p-n 
is directly calculable on the lattice, and more in- 
terestingly gD'D-w has been measured in CLEO 
through the D*+ -> D 0(+ \+^ decays || while 
B* — > Bit is prohibited kinematically. The rela- 
tion between g and gp'p^ is given by ChPT as 
gp'Piv = 2 j y/mp.mp/fn up to 0(l/M) and 
0{m%) corrections. Using an LSZ reduction of 
the pion and the PCAC relation, the calculation 
of above matrix element reduces to that of 



{Pip) I Ql^q I P*(p + q)), 



(26) 



namely, to the calculation of the form factors de- 
scribing the above matrix element. 

The first exploratory study to determine gp* p» 
on the lattice was made in the static limit by 
UKQCD |||. This year the determination of 
go'D-K at physical D meson mass was presented 
by Herdoiza in the quenched approximation pQj. 
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Figure 12. The 1/Mp dependence of gD-D-n and 
g in the chiral limit. 



> 

CL) 

O 









=~-i-~~i 


5 ' 


J l 


i 








l/a = 1 


.0 GeV I 



0.5 1 
|p| 2 (GeV) 2 



Figure [T^ shows the 1/Mp dependence of go* Dn 
and g in the chiral limit. It is found that the 
1/Mp dependence of the former is completely 
canceled by the factor yjmpmp. After interpo- 
lation to the physical D meson mass, they ob- 
tain gD*D-n = 18.8(2. 3)(^;J), which agrees with 
5/3*^=17.9(0.3) (1.9) measured by CLEO |§. 

It seems that the first step toward the determi- 
nation has been completed. The next step is the 
unqucnchcd simulation with special care to the 
chiral extrapolation. 

4. Improvements in methodology 

4.1. KS as light quark 

We have seen that the issue associated with chi- 
ral logarithm is present almost everywhere. One 
possible course to simulate light sea quark masses 
is to employ staggered fermions. 

At the last conference, Wingate et al. per- 
formed an exploratory study of B mesons with 
the NRQCD and staggered actions on a coarse 
lattice (a _1 =0.8 GeV), and pointed out a dif- 
ficulty associated with contaminations from un- 
physical modes [pl| . This year they repeated 
the similar calculations on slightly finer lattices 
using N f =0 (a" x =l GeV) and N f =3 (a' 1 =1.3 
GeV) configurations ^52|. They found that such 
unphysical modes can be essentially removed by 
the use of finer lattice. They also observed that 
the correlation functions with both parities is well 
separated by applying suitable fit forms and the 
constrained curve fitting |6j|. 



Figure 13. Comparison of the mass obtained from 
dispersion relation (open square) with that using 
perturbation theory (dashed line). The solid lines 
indicate the perturbative uncertainty in the lat- 
ter. 

They calculated B s meson spectrum, and made 

a comparison of Mp a obtained from dispersion 

relation with that using perturbation theory as 

shown in Fig. [13[ These two definitions of 

mass agree quite well. They also obtained 
^ Nf= °)=225±9(stat.)±20(pert.) MeV in thc 

quenched approximation and this value is reason- 
ably consistent with the current quenched world 
average /^ f=0) =200(20) MeV. More importantly 
they conclude that the unquenching dramatically 
improves the hyperfine splitting Mg* — Mb, and 
results in 42.5±3.7 MeV. These observations in- 
dicate that there is no practical problem in the 
combination of NRQCD and staggered fermions. 

4.2. Anisotropic lattice 

The motivation for the use of anisotropic lat- 
tices in heavy-light system is to obtain clean sig- 
nals in the form factor calculations. The method 
has been applied to the B — > 7r semileptonic 
decay by Shigemitsu et al. Q. The simu- 
lation is performed on a 12 3 x 48 anisotropic 
quenched lattice with a s /at=2.71 (l/a s =1.2GeV) 
using the tadpole-improved Symanzik gluon ac- 
tion, NRQCD for heavy quarks and D234 action 
for light quarks. They apply the constrained fits 
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Figure 14. Comparison of form factors obtained 
with iso- and aniso-tropic lattice given by Q. 

PI to the three-point functions to obtain the 
form factors shown in Fig. [l4|. Their results agree 
well with those with isotropic lattices. In spite of 
a relatively small number of configurations (199), 
the statistical error in this work is comparable or 
even smaller than those of JLQCD, which used 
NRQCD heavy quarks and was obtained with 
more than 1,000 configurations [j55| . 

4.3. Step scaling method 

At this conference a new method to calculate 
f B in non-perturbative accuracy was presented 
by Petronzio [psf , based on a non-perturbative 
recursive finite size technique. The explicit calcu- 
lation of f B is made in the Schrodinger functional 
setup with several quenched lattices. Most of rele- 
vant renormalization constants are obtained non- 
perturbatively. My understanding of this method 
is as follows. Let us start from the following iden- 
tity, 



/b(12 GeV)|i. 6 = / B (12 GeV)| . 4 

/s(12 GeV)| , 8 / B (12 GeV)|i. 6 
X / B (12 GeV)| .4 X /b(12 GeV)| . 8 : 



(27) 



where /b(12 GeV)|i denotes f B measured on the 
lattice with a _1 = 12 GcV and the physical volume 



L fm. The first factor /b(12 GeV)|o.4 is inter- 
preted as f B on a finite volume and the remaining 
factors are considered as corrections to obtain f B 
with large volume. It should be noted that with 
this lattice spacing the b quark can be suitably 
treated by the relativistic formalism and that the 
above equation holds even if the smallest lattice 
is in the deconfined phase. Since it is difficult to 
realize the lattice simulation with both a -1 =12 
GeV and L=1.6 fm simultaneously, consider the 
following, 



f B (12 GeV) 1 1.6 



f B (12 GeV)| .4 

f B (3 GeV)|i. 6 



f B (6 GcV)| . 8 
/b(6 GeV) I o.4 /s (3 GeV) | 



0.8 



(28) 



which is obtained by replacing the second and 
third factors by those obtained on coarser 
lattices but with the same physical volume. 
They performed above calculation and obtained 
/b=170(11)(5)(22) MeV and f Bs =192(9) (5) (24), 
where the first error is statistical and the second 
and third are systematic. These values agree well 
with the current quenched estimates (|^). 

The questions are how ( p8| ) is justified, what 
is f B in the deconfined phase, and what is the 
advantage of this method compared to the cal- 
culation on a single lattice with a -1 =3 GeV and 
1.6 fm. The possible systematic uncertainties are 
still unclear to me. The same calculation in the 
static limit would be helpful to make this clear. 

5. Summary 

We are now going into an exciting era because 
Belle, BaBar, Tevatron Run II and CLEO-c will 
be giving us a wealth of precision data for B and 
D mesons soon. 

One of the CKM matrix elements that these 
experiments should pin down is \V t d\- An impor- 
tant point realized recently is that the amplitude 
of B°-B° mixing necessary to convert the exper- 
imental mass difference to this matrix element 
suffers from a sizable uncertainty of 10-20% due 
to chiral logarithms in the chiral extrapolation. 
However, one can avoid this problem if one intro- 
duces the Grinstein ratio of decay constants. It 
is now possible to determine the relevant quanti- 
ties on the lattice as accurately as 5% and better, 
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and the precise determination of \Vtd\ is realized 
when AMb s , jo d and fo s are measured in the 
forthcoming experiments. 

For \V c b\, needed for the normalization of the 
matrix elements, the situation is more promis- 
ing because the form factor of the B — * D*lv 
decay can be determined on the lattice with a 
few percent accuracy. The lattice method for 
this calculation is established, and does not re- 
ceive significant effects from chiral logarithms. It 
only remains to apply the method to unquenched 
calculations, and independent checks using other 
processes are in progress. 

The extraction of \V u b\ from B — > n(p) is still 
theoretically challenging as it is so in experiment. 
The form factor calculations still have large un- 
certainties (~20%) even in the quenched approx- 
imation. Further improvements are necessary. In 
particular, it is essential to find out how one can 
secure clean signals from simulations and how to 
overcome the problem of chiral logarithms. 

To achieve better accuracy, improvements in 
the methodology are also crucial. For the is- 
sue of chiral logarithm, lowering the quark mass 
provides the direct route for resolution for which 
staggered light quarks may be helpful. Clean sig- 
nals in form factor calculations are obtained with 
anisotropic lattices. Finally the non-perturbative 
accuracy might be achieved by the development 
of the step scaling method. 
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